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Abstract 

Results for quantum mechanical calculations of the integral cross sections and corresponding ther- 
mal rate coefficients for para-/ortho-H2+HD collisions are presented. Because of significant astro- 
physical interest in regard to the cooling of primodial gas the low temperature limit of para- / ortho- 
H2+HD is investigated. Sharp resonances in the rotational state-resolved cross sections have been 
calculated at low energies. These resonances are important and significantly contribute to the 
corresponding rotational state-resolved thermal rate coefficients, particularly at low temperatures, 
that is less than T ~ 100K. Additionally in this work, the cross sections for the elastic HD+HD 
collision have also been calculated. We obtained quite satisfactory agreement with the results of 
other theoretical works and experiments. 
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I. INTRODUCTION 



We report rotational state-resolved cross sections and corresponding thermal rate coef- 
ficients for the fundamental elastic and inelastic collisions between hydrogen molecules in 
both ortho- and para- states and its isotopic partner, the HD molecule: 

ortho-/para-E 2 (j2) + HD(ji) -> ortho- /para-R 2 (j' 2 ) + HD(#). (1) 

Here H is hydrogen and D is deuterium. HD is the second most abundant primordial 
molecule after H2 and plays an important role in the chemistry of the interstellar medium 
fl, more specifically in rotational energy transfer collisions © fl fl Q. Also, HD itself is of 
interest as a light molecule with a small but constant dipole moment \d\ = (8.1 ± 0.5) ■ 10~ 4 
D |||. This ability of HD makes the study of this molecule an important undertaking in 
astrophysics and other fields, such as in condensed matter physics, see for instance [6]. 
Furthermore, Ho and HD molecules play a significant role in the formation of the primary 

un 

stars in the Universe [3, la], where the dynamics of H2+ H2, H2+ HD and HD+HD molecular 
collisions can make a substantial contribution to the formation of the primordial structure of 
stars |9|. The spectroscopic measurements from distant quasars with H 2 and HD molecules 
in the cold gases around them or on the line of sight can be very important and requires 
detailed knowledge of the cross sections at low energies 
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Although the field of molecular collisions between hydrogen molecules and its different 



isotopic combinations 



H 



research history 
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+Do, Do+HD, HD+HD and others, have already a rather long 
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there are substantial discrepancies between 
the results of different research groups. The recent appearance of the new potential energy 
surface for the H4 system makes it very useful to carry out new calculations for these 
important systems. 



By now the deuterium chemistry of the ear 



few research groups, see for example, [2J and 



Q and 



Jniverse has been extensively studied by a 



In these works the authors carried out 



quantum-mechanical state-resolved calculations for a wide range of temperatures. Schaefer's 
work |2| presents results for the thermal rate coefficients at very low temperatures down to 



~10 K. The corresponding cross sections have been calculated at collisional velocities ~ 12 
m/s, which are almost equal to ultra-low energies. 

The author of work [3]] used Schwenke's modified H2-H2 potential [ji], [3]. Schaefer 
applied a different surface, which was also a modification of older PESs, see for example 



2l| . However, as we already mentioned, the importance of the H 2 +HD rotational energy 



transfer collisions makes it very reasonable to carry out new calculations for the systems 



with the use of the BMKP P 
most accurate H 4 surface 22] 



IS published in 2002 which provides the latest and presumably 



In this work the scattering cross sections and their corresponding thermal rate coefficients 
are calculated using a non reactive quantum-mechanical close- coupling approach. In the 
work [4j] we already presented the details of the method we used, therefore in this work we 
continue beyond that point with the results for ort/io-H 2 +HD . The calculations are done for 
very low and low temperatures from 5 K to up to 300 K. Atomic units (e = m e — h — 1) 
are used throughout these calculations. 

II. RESULTS 

Before undertaking production calculations we carried out a large number of test calcu- 
lations to insure the convergence of the results with respect to all parameters that enter into 



the propagation of the Schrodinger equation 23]. The details of the quantum- mechanical 
computation approach are presented in 24j. The same test calculations were carried out 
in our previous works 25|, |26j for the ortho-/para-B. 2 +B.2 collisions, which involved the in- 
termolecular distance R3, the total angular momentum J of the four atomic system, the 
number of rotational levels to be included in the close coupling expansion Ni v i and others. 

In these calculations the VRTP mechanism has been used, which allows us to specify 
V(R\, R2, R3) explicitly rather than to expand the potential in angular functions, see the 



MOLSCAT manual [23]. As in our previous works we reached convergence for the integral 
cross sections, (j(j' 1 ,j 2 ;jij2,^), in all considered collisions. However, we have found in this 
work, that in the case of ort/io-H 2 +HD the convergence is slower than in the case of para- 
H 2 +HD. 

A. ortho-/para-~H.2-\-HD and HD+HD cross sections 

In Fig. 1 we present the results for the low energy ortho- / para-H.2 + HD and HD+HD 
elastic scattering collisions. In these calculations the collision energy ranges from 10~ 5 cm -1 
to 100 cm -1 . As can be seen at energies less than ~ 5 • 10 -4 cm -1 the elastic cross sections, 
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a e i, display almost constant values. This result is in good agreement with the non-relativistic 
quantum-mechanical scattering theory, that is at very low collision energies: 

Vei{Ecou ~ 0) = -nal = const, (2) 



where a® is the scattering length [27|]. From our calculations in the case of ort/io-H 2 +HD: 
o~ e i = tcclq = 380.70 -10" 16 cm 2 . The corresponding scattering length is ao{ortho-H.2 + HD) = 
11.01 ■ 10~ 8 cm. A rather sharp resonance is found in this cross section at energies around 
0.56 cm -1 . The value of this resonance is a res (ortho-B.2 + HD) = 1, 897.8 ■ 10~ 16 cm 2 . 

Below, in the same figure we also present the results from a para-H.2+HD scattering 
cross section. These results were not shown in Further, we computed the HD+HD 
elastic collision. These calculations are also done with the BMKP PES. Fig. 1 shows this 
result. It is seen, that the behaviour of all three cross sections is almost identical, although 
the values are rather different. Again, the collision energy ranges from 10~ 5 cm -1 to 100 
cm -1 . At energies less than ~ 5- 10~ 4 cm -1 the elastic cross sections display almost constant 
values: a e i(para-R 2 +HD) = 380.031- 10" 16 cm 2 , and cx e ; (HD+HD) = 274.84-10~ 16 cm 2 . The 
corresponding scattering lengths are: ao(para-H2+HD) = 11.0- 10~ 8 cm and ao(HD+HD) = 
9.35 • lO" 8 cm. 

Sharp resonances have been calculated at energies around 0.57 cm" 1 in the case of para- 
H2+HD and around 2.05 cm -1 for HD+HD. The values of these resonances are the following: 
a res {para-YL 2 + HD) = 1,888.84 ■ 10~ 16 cm 2 and a res (HD+HD) = 737.43 ■ 10" 16 cm 2 . 

Next, one can see from Fig. 1, that the elastic cross sections are closely related to each 
other and the calculated values of do's have reasonably similar values. However there are 
significant differences in the positions and the values of the resonances between ortho-fpara- 
H 2 +HD and HD+HD collisions. The reason for these differences is due to the permanent 
dipole moment of the HD molecules. This phenomenon is important at low energies and is 
responsible for the energy shift in the total Hamiltonian of the four-atomic systems. The 
additional interaction makes the shift of the position of the high value resonance in HD+HD 
elastic cross section relative to the positions of the resonances in the ortho-fpara-^+HD 
systems. 

In addition to the results of Fig. 1 in Fig. 2 our HD+HD elastic cross sections are shown 



together with the corresponding results from paper 



Please see Fig. 16 and Fig. 17 of 



13j . where theoretical and some experimental results for HD+HD are presented. It is seen, 



that the trend of the both cross sections is similar to each other, this is important, although 
the values are rather different at some collisional energies (velocities). We consider these 
results as a quite sufficient test of our calculations for the HD+HD collisions. 

Fig. 3 shows results for the state-resolved cross sections o~(j[j 2 ] jij 2 )(e) in the ortho-case 
of the collision (1). In the upper and middle plots we present our data for the (ji = 1, j 2 = 
f) — > (j[ = 0,j' 2 = 1) and (ji = 2, j 2 — 1) — > (j[ = j 2 = 1) quantum transitions together with 
;he results of Schaefer's work [2], which applied a modified version of Schwenke's potential 
2o| and a different dynamical quantum-mechanical method. In the lower plot results for 
(J i = 2, j 2 — 1) — > (j[ = 0,j 2 = 1) are presented together with the results of work [2]. The 
values of this channel cross section are lower by factor of 10 than in the case of the two 
previous cross sections. 

As can be seen for the transition quantum states (ji = l,j 2 — 1) — > (j[ = 0,j' 2 = 1) 
and (J i = 2,j 2 — 1) — > (j[ = j' 2 = 1) the agreement between our results and Schaefer's 
calculations is fairly good. However, for the (ji = 2, j 2 — 1) — > (j[ = 0,j' 2 = 1) channel cross 
section we found significant disagreements in the cross section. The difference between our 
cross section and Schaefer's results reaches the factor of 2 at a collisional velocity around 
270 m/s. 

In conclusion, in the case of ortho- / para-H. 2 -\-HD we reproduced all sharp resonances in 
the low energy region {4]. These high value cross sections should contribute significantly to 
the cooling process of the astrophysical media, therefore it is of crucial importance to carry 
out calculations for these cross sections accurately. However, in the case of higher transition 
states we found some disagreements (up to 100%) in the cross sections with work [2]. But 
these cross sections have values of relatively low magnitude. 



B. ortho- /par 0-H2+HD thermal rate coefficients 

In this section we present new results for thermal rate coefficients regarding the case 
of ortho- / par a-H. 2 +HD . Corresponding cross sections have been presented in the previous 
section and some of them have been published in work jj] . In we were able to reproduce 
all relevant resonances in the cross sections for different rotational channels j^, however 
significant differences and new resonances have been found in the cross sections with higher 
values of the rotational transition states. 



r; 



Fig. 4 shows temperature dependent thermal rate coefficients corresponding to two 
relevant cross sections from Fig. 2: for (ji = l,j 2 = 1) - ► (J[ = 0,^ = 1) and 
(ji = 2,j 2 = 1) — > j'2 — 1) quantum-mechanical rotational transition states. Again 

agreement with work |2| is very good. However there is significant disagreement with results 
of work 3]. 

Fig. 5 presents our results together with Schaefer's data from for the (ji = 1, j 2 = 
0) - U'l = 32 = 0), (ji = 2,j 2 = 0) - (j'i = 1, j 2 = 0) and & = 2, j 2 = 0) -> = j 2 = 0) 
thermal rate coefficients. One can see, that the first two thermal rates agree very well with 
the results of work [2], however the third (lower plot) is not in good agreement. As in Fig. 
3, there is significant disagreement with the corresponding results of work j^]. 



In J4] and in this work we provide cross sections and corresponding thermal rates for only 
relevant transition states, that is transitions with relatively high cross sections. Detailed 
calculations for many other transition quantum numbers and, probably, for collisions of 
other hydrogen isotopes should be done in the future. 

Through this analysis we can conclude, that the BMKP PES can reproduce the general 
behaviour of all considered cross sections and thermal rates in the ortho/para-B. 2 +HD col- 
lisions. For the lower quantum states we obtained sufficient agreement with the results of 
work [2!] and significant disagreement with [3]. Also, for transition states with higher values 
the BMKP PES provides rather small cross sections and thermal rate coefficients relative 
to the corresponding results of work 

III. CONCLUSION 

In this work the state-resolved close-coupling quantum-mechanical calculations for rota- 
tional excitation and de-excitation cross sections of the ort/io-/para-H 2 +HD and HD+HD 
collisions are presented. The latest, global BMKP surface for the H 2 -H 2 system has been 
appropriately adopted for H 2 +HD and HD+HD. The linear rigid rotor model for the H 2 and 
HD molecules is applied. A test of convergence and the results for cross sections with the 
BMKP PES are obtained for a wide range of kinetic velocities including very low values down 



to ~10 m/s. As in work [4] these results revealed, that the applied quantum-mechanical 
method together with the BMKP PES is able to provide reliable cross sections and thermal 
rate coefficients in collisions (1). 
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An interesting result of this work lies in the comparison between our calculations with 
the BMKP PES, calculations of j^], and with the results of 2]. As we mentioned, these 
works used pure quantum mechanical methods, but within the rigid rotor approximation: 
the distances between hydrogen atoms in the H 2 or HD molecules have been fixed. These 
works used different H4 potential energy surfaces. In Fig. 3 we showed that our results for 
the state-resolved cross sections &rd-^riv{Ecoii) are close to the corresponding cross sections 
of the relatively old work ['\. Furthermore, because the state-resolved thermal rate coeffi- 
cients k n i-+n'i'(T) are less sensitive to different collision parameters, for instance, interaction 
potentials, for k n i^ n 'y(T) we obtained even better agreement with the results of work 
This is seen in Figs. 4 and 5. 

However, the same type comparisons with the corresponding thermal rate coefficients 
from the newer work [3J demonstrated substantial disagreements, which are also seen in 



Figs. 4 and 5. The deflection between results of [3j and our data and Schaefer's results can 
reach almost 60%. In future works it will be very useful for comparison purposes to adopt 



other precise H 4 PES from the work 28j] and apply it to the ortho-/para-H 2 +HD collisions. 



The Deep and Johnson potential energy surface (DJ PES) was already successfully applied 



to H2+H2 in 2a, 12a, |22|. Also, it seems to us that it would be very useful to carry out 



new calculations for various differential cross sections da n i^ n ii>(6, ip)/dQ for the H2+H2/HD 



collisions as in 



30) with the use of different PESs. 
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Fig. 1 Elastic scattering cross section for ort/io-H 2 +HD, para-H 2 +HD and HD+HD at 
ultralow collsion energies. Calculations are done with the BMKP PES. 



Fig. 2 Elastic scattering cross section for HD+HD ca 
PES and the corresponding result from work 



culated with the use of the BMKP 

1 



Fig. 3 Rotational state resolved integral cross sections for ort/io-H 2 (j 2 ) + HD(ji) — > 
H 2 (j 2 )+HD(j|). Upper plot: initial states of HD and H 2 molecules are j\ = j 2 = 1 
and corresponding final states are j[ = and j 2 = 1. In the bottom plots: j\ = 2, 
j 2 = 1 and the corresponding final states are j[ — f 2 — 1 and j[ — 0, j 2 — 1. 
Calculations are done with the BMKP PES (bold lines), circles, diamonds and 

H 

triangles up are corresponding results from work |2J • 

Fig. 4 Rotational state resolved thermal rate coefficients for ortho-E. 2 (j 2 ) + HD^) — > 
H 2 (j 2 ) + HD(j^). In the upper plot the initial states of HD and H 2 molecules 
are ji = j 2 = 1 and the corresponding final states are j[ = 0,j' 2 = 1. In 
the bottom plot, ji = 2, j 2 = 1 and corresponding final states are j[ = j 2 = 1. 
Calculations are done with the BMKP PES (bold lines), triangles up and squares 
are corresponding results from work [2J]. Results from work [31] are also shown. 

Fig. 5 Rotational state resolved thermal rate coefficients for para-H 2 (j 2 ) + HD(ji) — > 
H 2 (j 2 ) + HD(jJ). Upper plot: initial states of HD and H 2 molecules are ji = 1 
and j 2 = and corresponding final states are j[ = j 2 = 0. In the bottom 
plots: ji = 2, j 2 = and the corresponding final states are j[ = 1, j 2 = and 



j[ = f 2 = 0. Calculations are done with the BMKP PES (bo 
triangles up and diamonds are corresponding results from work 
work Nl are also shown. 



d lines), circles, 
2|. Results from 
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FIG. 1: Elastic scattering cross section for ort/io-/para-ri2+HD and HD+HD at ultralow collsion 
energies. Calculations are done with the BMKP PES. 
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FIG. 2: Elastic scattering cross section for HD+HD calculated with the use of the BMKP PES 
and the corresponding result from work 
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FIG. 3: Rotational state resolved integral cross sections for ort/io-H2(j2) + HD(ji) — > H2U2) + 

HD(jJ). Upper plot: initial states of HD and H2 molecules are j\ = ji = 1 and corresponding final 

states are j[ = and j' 2 = 1. In the bottom plots: j\ = 2, j'2 = 1 and the corresponding final 

states are j[ = j' 2 = 1 and j[ = 0, j' 2 = 1- Calculations are done with the BMKP PES (bold lines), 

circles, diamonds and triangles up are corresponding results from work 
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FIG. 4: Rotational state resolved thermal rate coefficients for ortho-H.2^2) + HD(ji) — > H2(j 2 ) + 

HD(j^). In the upper plot the initial states of HD and H2 molecules are j\ = j% = 1 and the 

corresponding final states are j[ = 0,j' 2 = 1. In the bottom plot, ji = 2, ji = 1 and corresponding 

final states are j[ = j' 2 = 1. Calculations are done with the BMKP PES (bold lines), triangles up 

and squares are corresponding results from work [2(]. Results from work [3| are also shown. 
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FIG. 5: Rotational state resolved thermal rate coefficients for para-H2(j2) + HD(ji) — > H2U2) + 

HD(j^). Upper plot: initial states of HD and H2 molecules are j\ = 1 and j% = and corresponding 

final states are j[ = j' 2 = 0. In the bottom plots: j\ = 2, j'2 = and the corresponding final states 

are j[ = 1, j' 2 = and j[ = j' 2 = 0. Calculations are done with the BMKP PES (bold lines), circles, 

triangles up and diamonds are corresponding results from work Q]. Results from work 3] are also 
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